Tomonaga— Luttinger parameters and spin excitations 
in the dimerized extended Hubbard model 
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We study the one-dimensional extended Hubbard model with alternating size of the hopping 
integrals using the density-matrix renormalization group method. We calculate the spin gap, the 
Tomonaga-Luttinger parameter, and the charge-density-wave order parameter for various dimeriza- 
tions, interaction strengths, and band fillings. At half band-filling the spin and charge excitations 
are gapped but these gaps disappear for infinitesimal hole doping. At quarter filling, the Umklapp 
scattering in the half-filled lower Peierls band generates a gap for the charge excitations but the gap- 
less spin excitations can be described in terms of an effective antiferromagnetic Heisenberg model. 
Beyond a critical strength for the nearest-neighbor interaction, the dimerized extended Hubbard 
model at quarter filling develops a charge-density-wave ground state. The dimerization and the 
nearest-neighbor Coulomb interaction strongly reduce the Tomonaga-Luttinger parameter from its 
value for the bare Hubbard model. We discuss the relevance of our findings for the Bechgaard salts. 

PACS numbers: 71.10.Pm,71.10.Fd,78.30.Jw,72.15.Nj 



I. INTRODUCTION 

The Bechgaard salts are organic conductors which have 
attracted much interest over the last thirty yearsi^. 
Upon variation of the pressure, the temperature, and the 
anion X in (TMTSF) 2 X and (TMTTF) 2 X, these com- 
pounds exhibit a rich phase diagram, e.g., a supercon- 
ducting phase is found to lie in-between a paramagnetic 
metallic phase and a spin-density-wave phase. The sys- 
tems can be regarded as quasi one-dimensional due to the 
strong anisotropy of the transport along the three crys- 
talline axes. Recent experiments 3 support the view that 
the metallic phase can be characterized as a Tomonaga- 
Luttinger liquid at temperatures T > 100 K. Indeed, 
signatures of the Tomonaga-Luttinger liquid are the re- 
duced density of states at the Fermi energy as seen in 
angle-resolved photoemission spectroscopy^, the nega- 
tive temperature dependence of the c-axis resistivity 6 , 
the scaling behavior of high-energy range of the optical 
conductivity, the power-law temperature dependence 
of the Hall coefficient^, and the empirical relationship 
(TiT) -1 oc Xs(T) between the measured spin relaxation 
rate and the magnetic susceptibility in nuclear magnetic 
resonance measurements^ *^. Moreover, distinctly dif- 
ferent thermal conductivities for the charge and spin ex- 
citations have been reported which provide evidence for 
spin-charge separation^. 

All correlation functions in the Tomonaga Luttinger 
liquid display a power-law behavior with unusual, inter- 
action-dependent coefficients. Many of them are sim- 
ple functions of the so-called Tomonaga-Luttinger pa- 
rameter K p . Most experiments give K p « 0.2 for the 
Bechgaard salts. The single-band Hubbard model in 
which spin-1/2 electrons move on a chain and interact 
only locally is one of the best studied Hamiltonians for 
correlated lattice electrons. However, the model gives 
K p > 0.5 for all interaction strengths which shows that 



the long-range parts of the Coulomb interaction must be 
taken into account for a proper description of the Bech- 
gaard salts. In the extended Hubbard model the long- 
range parts of the Coulomb interaction are mimicked by 
a nearest-neighbor termi 3 ^*^. 

Other factors may also play an important role. For 
instance, the stacks of TMTTF and TMTSF molecules 
form dimerized chains and the alternation of the elec- 
tron transfer-matrix elements along the chain must be 
considered. Therefore, in this work we study the one- 
dimensional extended Hubbard model with alternating 
hopping amplitudes, i.e., the one-dimensional dimer- 
ized extended Hubbard model as the minimal one- 
dimensional, purely electronic model for the electronic 
excitations in the Bechgaard salts. The relevant bands 
in the TMTSF and TMTTF salts are filled with three 
electrons so that the system is quarter-filled in hole no- 
tation, and we use the hole picture in the following. 

There are few systematic studies of the dimerized ex- 
tended Hubbard model in the literature. Therefore, we 
investigate the model for various band fillings, with an 
emphasis on the vicinity of the commensurate fillings. In 
this way, our principle investigation of correlated elec- 
trons in quasi one-dimensional dimerized systems could 
be relevant also for other materials, e.g., for the inorganic 
spin-Peierls system CuGeO^i^. 

In our work we apply the density-matrix renormaliza- 
tion group (DMRG) method which is one of the most 
reliable numerical methods to study the low-energy prop- 
erties of one-dimensional correlated electron systems. 
Where applicable, we compare our results to the predic- 
tions from field theory and effective single-band Hubbard 
models. 

Our paper is organized as follows. In Sec. [0] we de- 
fine the dimerized extended Hubbard model and intro- 
duce the physical quantities of interest, namely, the spin 
gap, the charge-density-wave (CDW) order parameter, 
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and the Tomonaga-Luttinger parameter. In Sec. IIIII we 
separately present our DMRG results for the dimerized 
Hubbard model with and without the nearest-neighbor 
interaction, and discuss the experimental relevance of 
our investigations. We close with a short summary in 
Sect.lTvl 



II. MODEL AND METHOD 
A. Hamiltonian 



structure is shown in Fig. ^ The dispersion relation of 
the two Peierls bands is given by 



ei 2 (k) = ±Jt\+tl + 2tit 2 co$k for |fc| < — , (2) 
v 2a 

where a is the lattice spacing which we set to unity in 
the following. The gap between the two Peierls bands is 
A P = 2(t 1 -t 2 ). The total band width is W = 2(h+h). 
In the absence of a dimerization, for t\ — t 2 = t, we 
recover the band structure of the tight-binding model in 
the reduced zone scheme. 



In order to model the Bechgaard salts, we focus on 
the transport of a chain of stacked molecules and regard 
a single TMTTF or TMTSF molecule as a site. The 
chain has a geometrical (Peierls) modulation. Besides 
the intra-molecular Coulomb interaction, we should take 
into account a nearest-neighbor Coulomb repulsion be- 
cause of the fairly short inter-molecular distance. Thus, 
our model Hamiltonian of choice is the one-dimensional 
dimerized extended Hubbard model for spin- 1/2 elec- 
trons on L lattice sites 



Z,odd 



/.even 



V^2{ni-n)(h l+1 -n) , (1) 



where c, {cia) is the creation (annihilation) operator of 



an electron with spin a =T,i at site I, hi a — &i a c-io is 
the number operator, and hi = n/j + hi[. The total 
number of electrons is N = + N±, and n — N/L 
is the average number of electrons per lattice site. The 
electron transfer matrix elements t\ and t 2 < t\ model 
the dimerization of the chain, U is the strength of the 
Hubbard interaction, and V parametrizes the nearest- 
neighbor Coulomb repulsion. We call a pair of sites which 
is connected by the hoping amplitude t\ a 'dimer'. 
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FIG. 1: Band structure for non-interacting electrons in a 
dimerized chain. 

The dimerization splits the tight-binding cosine band 
into a bonding band ('lower Peierls band') and an anti- 
bonding band ('upper Peierls band'). The bare band 



B. Physical quantities 

In this work we employ the DMRG method which pro- 
vides very accurate data for ground-state properties of 
one-dimen sional correlated electron systems; for a re- 
view, see |!7ll8l . We use the DMRG to calculate the 
spin gap A s , the CDW order parameter Xi an d t ne 
Tomonaga-Luttinger parameter K p . To this end, we 
consider a chain with L/2 dimers with L/2 even for 
a two-band system. We study chains with up to 320 
sites and open-end boundary conditions. We keep up to 
to = 3600 density-matrix eigenstates in the DMRG pro- 
cedure and extrapolate the calculated quantities to the 
limit m — > oo. In this way, the maximum error in the 
ground-state energy is below 10~ 6 ii. Lastly, we extrap- 
olate our finite-size results to the thermodynamic limit, 
L — > oo. 

The spin gap is defined by 

A s = lim A S {L) , 

L >oc 

A S (L) = ^(L,J\T T + l,.Nj-l)-^(£,tfT,W;)>(3) 

where Eq(L, iV-|-,iVjJ is the ground state energy of a sys- 
tem of length L with JV-f up-spin and N± down-spin elec- 
trons. 

Later in this work, we shall focus on the CDW ground 
state of our model JIJ at quarter band filling. For large 
enough nearest-neighbor repulsion V we expect a CDW 
with a wave vector Qcdw = 4fcp. Here, fcp = 7rn/2 is the 
Fermi wave number. At quarter band filling, n = 1/2, we 
have kp — ir/A which corresponds to a half- filled lower 
Peierls band. 

The order parameter for the 4fcp-CDW phase is defined 

by 



X = lim x{L) , 

L—^oo 



X(L) = 



r + 2 



(L+r)/2+l 

E (-!)'<% 

l=(i-r)/2 



(4) 



(5) 



In |J5J the summation over the lattice sites I is restricted 
to a region r around the central site of the chain in order 
to reduce the edge effects. We set r = 2 for a systematic 
extrapolation to the thermodynamic limit. Of course, 
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the extrapolated results should be independent of the 
choice of the range r. On finite lattices and for open- 
end boundary conditions, the Friedel oscillations from 
the edges result in a finite value for x(L), and a well- 
controlled finite-size extrapolation is mandatory. 

For the calculation of the Tomonaga-Luttinger param- 
eter K p we use a new method which we proposed re- 
cently—. The Tomonaga-Luttinger parameter K p deter- 
mines the long-range decay of the density-density correla- 
tion function in the metallic Tomonaga-Luttinger liquid 
ground state. It is defined by the ground-state expecta- 
tion value 



1 L 

C NN (r) = ~£)<n J+r n,) - (n 1+r )(n,) 



(6) 



i=i 



Using conformal field theory it can be shown 20 - 21 that 
the asymptotic behavior for 1 <C r -C L is given by 



C NN (r) ~ - 



K p A cos(2k F r) 



7rr) 2 



ln~ 3/2 (r) + 



(7) 



22.23,24.25 



where A is a constant. In previous approaches 
K p was extracted from the Fourier transformation 
of C NN (r) but in a real-space DMRG approach the ac- 
curacy of the correlation function becomes increasingly 
worse as the distance r increases. In Ref. we calcu- 
lated the density-density correlation function directly in 
Fourier space. We address 



N(q) 

where n(q) is given by 



2 
L 



n(q)n(-q)) 



n(q) =Y,e l 



i(g/2)(J+l/2-r c ) 



(8) 



(9) 



2. odd 



Here, r c = (L + l)/2 denotes the central position of the 
chain. The derivative of N(q) at q = directly gives the 
Tomonaga-Luttinger parameter. In practice, we obtain 
it from 



K, 



p = lim Kp(L) 



K P {L) 



L- 

L 



-N 



(10) 



For a precise calculation of K p is important to tar- 
get not only the ground state |$o) but also the state 
= n(-q)\9 Q ) in the DMRG procedure; see Ref. [Hj 
for further details. 

The Tomonaga-Luttinger parameter is well defined 
only for the metallic Tomonaga-Luttinger liquid. Later 
we shall investigate K p for insulators which are infinites- 
imally doped away from their commensurate doping n c . 
In these cases we give 



K p (n — > nf) = lim K p in = n c ± — 



(11) 



This approach is very successful for the single-band Hub- 
bard model, as demonstrated in Ref. 0|. 



C. Effective models 

For not too small dimerizations, ^Ai ^ 0.9, we can 
map the dimerized extended Hubbard model to an effec- 
tive single-band extended Hubbard models. The upper 
Peierls band can be integrated out and we are left with 
a Hubbard chain with L/2 dimer sites Id with effective 
parameters, 

H cS = ioff^( a L+i CT ai ^+ n - c -) + c/ cffX!^dT^dl 

Id ^d 

+V eS J2(ni d -l)(n u+1 -l), (12) 
id 

= ( ia ) 



^^.^FSMzM, (M) 



v cS 



V 



(15) 



The band filling is n e ff = 2n so that /cF.cff = t^n and 
^F.eff = i2sin(7rn). Note that U c s/t e s can be large even 
when U/ti is small, e.g., U e s/t g = 8.8 for U = ti, V = 0, 
and t 2 /h = 0.1. 

For V < V c the quarter-filled dimerized extended Hub- 
bard model describes a Mott-Hubbard insulator with 
gap-less spin excitations. Inn this parameter region, the 
spin degrees of freedom of the effective single-band Hub- 
bard model i|12fl can be described by an effective Heisen- 
berg model, 



-ffhcis.cff — </cff^^Si d ■ S/ d + i, 



(16) 



where S; d is the spin operator for a dimer located at 
position Id- Up to second-order in t2/U e s, we have 



J e s(V) = 



4*2 



8*i + 2U + V - 2y/(U- V) 2 + m'( 



(17) 



III. RESULTS 

A. Dimerized Hubbard model 

First, we consider the dimerized Hubbard model, i.e., 
we set V = in Q . 



1. Tomonaga-Luttinger parameter 

In order to demonstrate the accuracy of our method, 
we address the Tomonaga-Luttinger parameter at small 
interactions, U < W, as a function of the dimerization 
in the metallic regime, n = 0.4. To lowest order in the 
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FIG. 2: Tomonaga-Luttinger parameter K p from the DMRG 
approach (symbols) in comparison with the predictions from 
the g-ology method (solid lines), as a function of U/ti for 
t 2 /ti = 1, 0.9, 0.5, 0.3, 0.1 (from top to bottom) at n = 0.4 for 
the dimerized Hubbard model. 



couplings, g x = g 2 = 9a = U/2, 
ology' approach predict a 2 ^ 2 ? 



the field-theoretical 'g- 



2itvf 



2ttv f + U 



where the Fermi velocity v-p is given by 



v F 



t\ti sin fep 



Vff+lf+^i^ cos(fc F ) 



(18) 



(19) 



This result can be systematically improved with the func- 
tional Renormalization Group method 2 ^. 

In Fig. |21 we compare the Tomonaga-Luttinger param- 
eter as calculated from the DMRG approach, eq. (|10fl . 
to the g-ology prediction IjlH)). We plot K p as a func- 
tion of U/ti for various dimerization strengths t 2 /t\ at 
band filling n — 0.4. The system is metallic for all in- 
teraction strengths. For all dimerizations, K p decreases 
monotonically with increasing Coulomb interaction and 
finally approaches K p (U — > oo) = 1/2, as expected from 
the nondimerized Hubbard model. For small dimeriza- 
tion, t 2 /ti > 0.5, the DMRG results agree very well with 
those from the g-ology approach for all U < W . For small 
U/ti, K p decreases weakly and monotonically with t 2 /t\. 
This can be understood from the corresponding decrease 
of the bandwidth, W = 2(ti + t 2 ), with a corresponding 
reduction of the Fermi velocity. 

When the dimerization is large, £2/^1 ^ 0.5, and the 
Hubbard interaction is large, U > W/2, the results from 
g-ology substantially deviate from the numerically exact 
DMRG results. The Tomonaga-Luttinger parameter K p 
decreases rapidly with decreasing t 2 /ti, and the g-ology 
predictions quickly violate the constraint K p > 1/2. Ap- 
parently, higher-order corrections in U /W beyond the 
one-loop calculations needed to be considered. 

As our next application, we investigate the Tomonaga- 
Luttinger parameter as a function of the band filling and 
the interaction strength. In Fig. [3 we show K p from the 




FIG. 3: Tomonaga-Luttinger parameter K p as a function of 
the band-filling n for various dimerizations: (a) ti/ti = 0.9, 
(b) <2/ti = 0.5, and (c) tijtx = 0.1. In each figure, U/ti = 
1,2,6 from top to bottom. Open circles denote the DMRG 
results in the dimerized Hubbard model, and dotted lines are 
guides for eyes. Solid lines give the exact result for the single- 
band Hubbard model with hopping integral t = (ii + i2)/2. 



DMRG method as a function of n for various interac- 
tion strengths U/ti and dimerizations: (a) tz/ti = 0.9, 
(b) t2/ti = 0.5, and (c) t 2 /ti = 0.1. For comparison we 
also plot the exact results for K p from the Bethe Ansatz 
for the one-dimensional single-band Hubbard model with 
the same band width, t = (t± + t 2 )/2. 

When the dimerization is small, £2/^1 = 0.9, we again 
find a good general agreement between the results for the 
dimerized Hubbard model and the single-band Hubbard 
model with the same total bandwidth. An exception is 
the narrow range around quarter band filling, n = 1/2. 
At quarter filling, the lower Peierls band is half filled and 
the Umklapp scattering becomes a (marginally) relevant 
perturbation which turns the metallic phase into a Mott- 
Hubbard insulator where K p is not well defined, and we 
give the value for infinitesimal doping, see eq. I|llf> . As 
expected from field theorySL 2 ^, and confirmed numeri- 
cally, we have 



1± 

2 



(20) 



for the density-driven Mott transition for all interaction 
strengths. This follows from the mapping of the quarter- 
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filled dimerized Hubbard model to the effective single- 
band Hubbard model at half band-filling. Therefore, K p 
strongly changes as a function of density in the vicin- 
ity of quarter filling even for small dimerizations. The 
effect becomes more prominent with increasing dimer- 
ization strengths, see Fig. 

When the dimerization is large, t^/h = 0.1, the single- 
band Hubbard model does not provide a good start- 
ing point for the analysis anymore. Instead, for large 
ti/t 2 we rather consider the effective single-band Hub- 
bard model (|12fl for V = 0. Because of the strong effec- 
tive on-site interaction U e g/t e s, the Umklapp scattering 
strength becomes very large. For instance, the effective 
couplings at t 2 /ti = 0.1 are estimated from eq. lfTI|l as 
U eS /t eS = 8.8,15.3,27.9 for V = and U/h = 1,2,6, 
respectively. Therefore, the values for K p are rather 
small for all U /t\ > 1. Moreover, the effective single- 
band Hubbard model always gives the correct result 
K p (n — l/2 ± ) = 1/2 because the quarter-filled dimer- 
ized Hubbard model maps onto the half-filled single-band 
Hubbard model which describes a Mott-Hubbard insula- 
tor for all interaction strengths. 

As seen in Fig. 21 the quantitative agreement for K p 
from the dimerized Hubbard model and the effective 
single-band Hubbard model is quite good for all U/h at 
t2/ti = 0.1. Note that the effective Hubbard model dis- 
plays its particle-hole symmetry around n = 1/2 which 
the dimerized Hubbard model obeys only for t 2 /h — * 
or U jt\ — > 00. 

0.7 



T 



T 

U/h = 1 



T 



T" 



— U e g/t eS = 8.8 

U/h = 2 U eS /t eS = 15.3 

□ U/h = 6 U eB /t eS = 27.9 




FIG. 4: Tomonaga-Luttinger parameter K p from the DMRG 
approach for large dimerization, t^/tx — 0.1, in comparison 
with the analytical result for the effective single-band Hub- 
bard model. Recall that the band filling n of the dimerized 
Hubbard model corresponds to a filling 2n for the effective 
single-band Hubbard model. 



2. Spin excitations 

As our second quantity of interest we study the spin 
degrees of freedom at and around some commensurate 
band fillings. At half filling, n — 1, the dimerized Hub- 
bard model is a band-Mott insulator for all U /t\ > 0, and 
we expect a finite gap for spin excitations for all U /h . 



For small interaction strengths, the spin gap is of the or- 
der of the Peierls gap, A s (U/h ->• 0) = A P = 2(h - i 2 )- 
For large interactions, the spin degrees of freedom of the 
dimerized Hubbard model can be described by the one- 
dimensional Peierls-Hubbard model so that the spin gap 
to lowest order in ti/U becomes 



-it'i 



ti-tf 



2/3 



(21) 



in accordance with the results for the corresponding 
Peierls-Heisenberg model^S.. Eq. I|21(l is applicable for 
U/h > 4. In the inset of Fig. [5] we show two ex- 
amples for the finite-size scaling of the spin gap J3J), 
(t 2 /h = 0.5, U/h = 10) and (t 2 /h = 0.9, U/h = 5). 
The dependence of the gap on the system size is quite 
small because in the ground state individual spin sin- 
glets are formed on the dimers so that the gapped spin 
excitations are rather localized in space. 



0.02 



0.01 




0.015 



FIG. 5: Extrapolation of the spin gap A S (L) of the dimer- 
ized Hubbard model. Solid symbols represent the results for 
weak dimerization (ti/t\ = 0.9, U/ti — 5) at infinitesimal 
doping of the band-Mott insulator (n — 1, triangles) and at 
small doping (n — 0.95, circles). Open symbols give the re- 
sults for intermediate dimerization (t-2/ti — 0.5, U/ti = 10) 
at infinitesimal doping of the band-Mott insulator (n = 1, 
triangles), at the electron densities n = 0.95 (circles), n — 0.8 
(squares), and at infinitesimal doping of the Mott-Hubbard 
insulator at quarter filling (n = 0.5, lower triangles). 
Inset: Extrapolation of the spin gap of the band-Mott insu- 
lator at half band-filling for (ti/t\ = 0.9, U/ti = 5) (solid 
diamonds) and for {t-2./t\ — 0.5, U/ti = 10) (open diamonds). 



It is more interesting to study the doping dependence 
of the spin gap. In Fig. [5] we plot A S (L) as a func- 
tion of system size for (t 2 /h = 0.5, U/h = 10) and for 
(h/h = 0.9, U/ti — 5) for several band fillings. As seen 
from the figure, the spin gap vanishes for all electron den- 
sities. In particular, at half band-filling it disappears as 
soon as the system is doped with an infinitesimal amount 
of holes. This can be understood in terms of the spin 
excitations of a half-filled system with two holes. Let 
us assume that the two holes are confined to a dimer. 
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Then, a spin excitations would remain the same local ex- 
citation as in the perfectly half-filled system which costs 
the finite energy l|21|l. However, the holes are actually 
delocalized over the system because the breaking of two 
spin dimcrs cost twice A s but the gain in kinetic energy 
is approximately 



E it ~ 2(ti - t 2 ) 



(22) 



E it is always larger than 2A S . The mobile holes leave 
behind at least two broken spin dimers whose spin exci- 
tation energy vanishes in the thermodynamic limit. 

Apparently, the dimerized Hubbard behaves differently 
from the two-leg Hubbard ladder at half band-filling 
where a spin-singlet pair is formed on each rung. The 
spin gap in the ladder system remains finite for finite 
hole doping. There the spin-singlet pairs themselves are 
mobile so that in the ground state an additional pair of 
holes is actually confined to a rung because the gain in 
kinetic energy due to the hole motion is smaller than 
the combined loss in the pairing energy and the kinetic 
energy of the spin dimcrs. 

Finally, we investigate the spin gap for the quarter- 
filled dimerized Hubbard model at infinitesimal doping. 
In Fig. El we plot the size-dependence of the spin gap 
for the infinitesimally doped Mott-Hubbard insulator at 
quarter filling for (t^/ti = 0.5, U/t\ = 10). The extrapo- 
lated values are zero for all dimerization and interaction 
strengths. Therefore, the spin-gap liquid, suggested in 
the one-dimensional dimerized t-J model- 1 is not real- 
ized in the dimerized Hubbard model. 



B. Dimerized extended Hubbard model 

Now we turn to the case V ^ in Q). We focus on the 
region around quarter filling where the nearest-neighbor 
interaction can lead to a CDW phase. This is known for 
the extended Hubbard model whose ground-state phase 
diagram was studied in detail recently^i^. 



1. Charge order 

Previous studiea22i2iL2i suggested that the presence of 
a dimerization suppresses the CDW phase. Therefore, 
we investigate the dependence of the critical coupling V c 
for the onset of the CDW. To this end we calculate the 
CDW order parameter \ from J5| as a function of the 
dimerization strength. For V = we have \ = whereas, 
for large V, the CDW order parameter approaches its 
classical value, x{V — * °°) = 0-J5- 

In Fig. we show the order parameter x(V) as a func- 
tion of V/ti for dimerizations i 2 Ai = 1,0.9,0.5,0.1 for 
fixed U/ti = 10 at quarter band filling. In the ab- 
sence of a Peierls modulation, t^/ti = 1, i-e., in the ex- 
tended single-band Hubbard model, x^Y) 1S finite above 
V c /ti ks 2.65, in agreement with previous work^i*!*^. 



X 



0.5 
0.4 
0.3 
0.2h 
0.1 



~~ 1 r 
□ hlh 

x i 2 Ai 




0^-s-o-43H i k m 

2 3 



FIG. 6: CDW order parameter x extrapolated to the ther- 
modynamic limit L — > oo for t^/ti = 1, 0.9, 0.5, and 0.1 with 
fixing U/ti — 10 at quarter filling. Lines are guides to the 
eyes. Inset: Estimated critical interaction strength V c /ti for 
the CDW transition as a function of tijt\. 



Apparently, the dimerization enhances the charge fluc- 
tuation on each dimer, and, consequently the tendency 
towards charge order is reduced. 

In the presence of a dimerization the critical value for 
the onset of the CDW increases with increasing dimer- 
ization. Moreover, x(V> W*i < 1) r i ses U P sharply above 
KfeAi) even when t 2 /ti is close to unity. We speculate 
that the transition remains continuous for all finite i 2 Ai 
but the slope is infinite for all i 2 Ai > 0- In the inset of 
Fig.[S]we show the critical value V c /ti as a function of the 
dimerization strength ia/ii- We find that V c /t\ changes 
rapidly for small i 2 Ai an d quickly saturates at its classi- 
cal value for i 2 Ai = 0. The value V c (t 2 /ti = 0) = 4ti is 
readily explained by considering an isolated dimer. In the 
isolated-dimer limit the energies of the Mott-Hubbard 
insulator and the CDW insulator are 



Ef u /L 

~/L 



pCDW 







-h + V cS = -h + V/A , (23) 

(24) 



so that the criterion for the (discontinuous) transition is 
£ MH (V C ) = E^ W {V C ) which immediately gives V c /h = 
4. 



2. Tomonaga-Luttinger parameter 

In the absence of a dimerization, the Tomonaga-Lut- 
tinger parameter decreases as a function of V/ 1 for fixed 
U/t > 4 and reaches K p — 0.25 at the critical coupling. 
When the CDW insulator is infinitesimally doped the 
system metalizes and K° DW = l/8 19 i 29 i 36 . 

For a finite dimerization, the quarter-filled system is a 
Mott-Hubbard insulator for small V/ti and finite U/t\. 
At infinitesimal doping we find K™ H (V < V c ) = 1/2 
below the transition, independent of V . This is readily 
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understood from the fact that the effective model is the 
extended single-band Hubbard model at half band-filling 
for which the field-theoretical arguments for a density- 
driven Mott transition still apply. A qualitatively and 
quantitatively different behavior emerges from the tran- 
sition to the CDW insulator at V c . The Tomonaga- 
Luttinger parameter drops from K p = 1/2 in the in- 
finitesimally doped Mott-Hubbard insulator to K p < 
1/8, as we shall discuss in more detail now. 

The dimerization has two prominent effects on K p . 
First, it increases the strength of the Umklapp scatter- 
ing which makes K p smaller. Second, the dimerization 
suppresses the CDW instability which tends to make K p 
larger. These effects are most apparent around quarter 
filling where the two tendencies compete with each other 
close to the CDW instability. Both effects increase upon 
decreasing t^/t-^. The first effect continues to develop 
progressively and leads to U c s/t c e — > oo as t%jtx — * 0. 
As shown in Sect. Iff! B II the second effect develops fast 
as a function of the dimerization and quickly saturates. 
Therefore, we expect that the first effect, a reduction of 
K p upon dimerization, is more prominent but for quarter 
filling and in the vicinity of the transition to the CDW 
phase. 

The reduction of K p with dimerization can actually be 
inferred from the g-ology approach where the Tomonaga- 
Luttinger parameter near quarter filling is given by 




12 3 4 12 



FIG. 7: Tomonaga-Luttinger parameter K p in the dimer- 
ized extended Hubbard model as a function of the nearest- 
neighbor Coulomb interaction V/ti for U/ti = 6 and various 
dimerizations: iz/ti = 1 (dashed line), <2/ti = 0.9 (filled 
circles), £2/^1 = 0.7 (open triangles), £2/^1 = 0.5 (crosses), 
t-zjt\ = 0.3 (filled triangles), and t<i/t\ = 0.1 (open squares). 
The band filling is (a) n = 0.475 and (b) n = 0.525. 



K n 



2ttu f - V 
2nv F + U + 5V 



(25) 



with vf — txt^j \/t\ + 1\. The formula shows that K p 
decreases monotonously as a function of V and oit 2 /tx- 
Naturally, g-ology cannot cover large dimerizations or 
the transition region where the increase of K p (V) upon 
dimerization becomes apparent. 

In Fig. we present the DMRG results for K p as a 
function of V/ti at a hole doping and an electron dop- 
ing of 5%, (n) = 0.5 ± 0.025, for U/h = 6 and vari- 
ous dimerizations. The numerically exact DMRG results 
confirm the general expectations as expressed in the g- 
ology formula 1)25(1. The Tomonaga-Luttinger parame- 
ter decreases monotonously with V/tx for all dimeriza- 
tions and, in general, it decreases as a function of t 2 /tx 
for fixed V/t\. The fact that K p is almost independent 
of t%/tx for fixed 1 < V/tx < 2 can be attributed to 
the above-mentioned competition between the Umklapp 
scattering and the charge ordering. For certain parame- 
ter regions, a change in the dimerization strength has al- 
most no net effect on K p because a change in the strength 
of the Umklapp scattering is compensated by a change in 
the mobility of the charge carriers. For the same parame- 
ter set (V/t x ,t 2 /tx), K p is generally somewhat smaller for 
the hole-doped case than for the electron-doped case but 
there is no difference in the qualitative behavior. This 
had to be expected because the system is particle-hole 
symmetric around quarter filling to lowest order in t 2 /t X - 
From now on we shall focus on the case of hole doping. 



The Tomonaga-Luttinger parameter K p (V) changes 
most rapidly in the region 2 < V/tx < 4 where the 
quarter-filled system undergoes the charge-ordering tran- 
sition. For V/tx > 4, we can interpret the system as a 
doped CDW insulator. In this region, we find that the de- 
pendence of K p on the nearest-neighbor interaction V/tx 
is much weaker. This can be understood from the Tay- 
lor expansion of K p for a slightly doped CDW insulator. 
Above the transition point (V > V c ) we generally ex- 
pect^ that for 6 — 1/2 — n <C 1 we have 



K p (t 2 , U, V,l/2-S) = Kf DW (t 2 , V) + 



h(t 2 ,U,V) 



(26) 

where t\ is used as energy unit. The prefactor h(t 2 , U, V) 
diverges exponentially at the critical interaction strength 
V c but it rapidly tends to a constant for large V. 

For infinitesimal doping, the Tomonaga-Luttinger pa- 
rameter of the CDW insulator i^ DW (i 2 /£i, V/t x ) also 
displays a smooth behavior as a function of V/tx and 
t 2 /h. In Fig. El we show K° DW (t 2 /tx,V/tx) for U = 00 
and various dimerizations. As for the case of a finite 
doping we see that the dimerization tends to reduce 
the Tomonaga-Luttinger parameter. In the CDW phase 
this tendency is somewhat compensated by the influence 
of the nearest-neighbor Coulomb interaction which, for 
large interactions and for small doping of the CDW state, 
delocalizes the holes over the system and therefore in- 
creases the charge fluctuations which determine K p via 
eq. 0. 
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relations are dominant for slightly doped Mott-Hubbard 
and CDW insulators, respectively, and their correlation 
functions decay algebraically with the asymptotical be- 
havior C* 2fcF - SDW (r) ~ r - 1 - K " for Mott-Hubbard insu- 
lators and C 4feF - CDW (Y) - r~ 4K " for CDW insulators. 
Thus, the value K p = 1/3 discriminates the two types of 
(doped) insulators at finite doping. For an infinitesimally 
doped insulator we correctly find K p (V = Vd c ) = 1/3 for 
Vdc = V c but Vdc becomes actually smaller upon doping, 
as seen in Fig. |5J 



FIG. 8: Tomonaga-Luttinger parameter for the infinitesi- 
mally doped CDW insulator K° DW as a function of V/t\ for 
t 2 /ti = 0.9, t 2 /ti = 0.5, and t 2 /*i =0.1 at U/h = oo. The 
solid line corresponds to K^ I3W = 1/8 when the dimerization 
is absent (t 1 = t%), and the dotted lines are guides for the 
eyes. 

The most important observation is the magnitude of 
the Tomonaga-Luttinger parameter for the doped insu- 
lators. For infinitesimal doping we find K p — 1/8 in the 
absence of dimerization and even K p < 1/8 in the pres- 
ence of a dimerization. These small numbers persist for 
finite doping, as seen in Fig. \7\ Therefore, depending on 
the choice of the dimerization and the nearest-neighbor 
Coulomb interaction, one can easily find parameter re- 
gions where 0.1 < K p < 0.3 can be realized for slightly 
doped quarter-filled chains. 



~ a few % hole doping 5 ~ a ^ ew ^ electron doping 




V 



FIG. 9: Schematic phase diagram of the dimerized extended 
Hubbard model around quarter filling as a function of n and 
V. Variations of K p are displayed by contour lines. Darker 
(brighter) color denotes smaller (larger) values of K p . 

Fig. |2| summarizes our findings for the Tomonaga- 
Luttinger parameter in a schematic phase diagram for the 
slightly doped quarter-filled dimerized extended Hubbard 
model. The Mott-Hubbard insulator (CDW insulator) 
can be characterized as 2fcp-SDW (4rCp-CDW) states at 
quarter filling. Hence, the 2rCp-SDW and 4/cf-CDW cor- 



C. Comparison with experiment 

At last, we compare our theoretical result with experi- 
ments on (TMTTF)2X. The electron transfer matrix ele- 
ments^ are estimated to be (ti,^) = (137meV,93meV) 
for X=PF 6 , (h,t 2 ) = (140meV, lOOmeV) for X=C10 4 , 
and (ti,t 2 ) = (133 meV, 119meV) for X=Br, i.e., t 2 /h = 
0.68,0.71,0.89 for X=PF 6 , C10 4 , Br, respectively. From 
the comparison with the optical gap 3 -' 39 the Coulomb 
parameters are estimated to be U/t\ ~ 7.0 and V/ti ~ 
2.8 for (TMTTF) 2 PF 6 . A comparison with Fig. □ 
and [5] shows that this parameter set leads to K p w 
0.25, in agreement with experimental estimates for the 
Tomonaga-Luttinger parameter from the temperature 
dependence of the resistivity 6-40 . In view of the CDW 
state observed below T « 100 K 41 i 42 i 43 , the nearest- 
neighbor interaction could be even larger than V jt\ = 2.8 
which would further reduce K p . 

Unfortunately, such values for the nearest-neighbor 
interaction V/t± appear to contradict the results for 
the effective exchange interaction as deduced from the 
high-temperature data from the electron-spin-resonance 
(ESR) measurements^, J exp = 420 K, 430 K, 500 K for 
the anions X=PFg, CIO4, Br, respectively. In the pres- 
ence of the dimerization and at quarter band-filling we 
can start from the effective extended single-band Hub- 
bard model (|12fl and the spin degrees of freedom can be 
described in terms of the effective Heisenberg Hamilto- 
nian (|16fl . For U jt\ = 7.0, the bare Hubbard model, 
V = in U3), gives J cS (V = 0) = 499 K, 564 K, 841 K. 
The good agreement of the experimental and theoretical 
data for V = implies that the nearest-neighbor inter- 
action ought to be rather small. In particular, the a 
value V = 2.8ti for (TMTTF) 2 PF 6 , leads to J cS (V = 
2.8ii) = 222 K, a factor of two smaller than the exper- 
imental estimate. Additinally, with small V to adjust 
^oxp, the resulting theoretical prediction for K p w 0.5 
from Fig. is not compatible with the experimental es- 
timate, 0.2 < K p < 0.3. 

In order to reconcile this discrepancy we note that, in 
the ESR measurements, the curves are fitted to provide 
a good agreement with the Eggert-Affleck-Takahashi 
model 4 *! for the spin susceptibility of the S — 1/2 an- 
tiferromagnetic Heisenberg chain at elevated tempera- 
tures. However, substantial deviations occur for small 
temperatures, T < 100 K. They could be the result 
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of a dimensional crossover— and the transition to the 
CDW phase. We are tempted to attribute the devia- 
tions to an effectively larger nearest-neighbor interaction 
at low temperatures. Recall that our electronic model is 
purely one-dimensional, and neither covers the influence 
of phonons^£*^ nor does it give an account on the screen- 
ing of the electron-electron interaction which may change 
drastically in the vicinity of the transition to the CDW 
state. Therefore, temperature may have a quite substan- 
tial influence on the value of the effective ^-parameter 
in our model so that eq. Ijl7|l cannot be applied with the 
values for V/t\ at T = to explain the susceptibility data 
forT>100K. 

In fact, in the CDW phase, the effective exchange 
interaction is given by J ^ DW jt x « 4i|/(2[/y 2 ) which 
results in J^jp w — 14 K if we use the parameters for 
(TMTTF) 2 PF 6 . If the spin susceptibility could be mea- 
sured in the (one-dimensional) CDW phase, the exchange 
interaction should be an order of magnitude smaller than 
in the high-temperature phase. 

In (TMTSF) 2 PF 6 , the hopping amplitudes are esti- 
mated as (<i = 252 , meV, t% = 209 mcV) and the effec- 
tive Coulomb interactions are found to be weaker, U /t± ~ 
5. Again, a weak nearest-neighbor Coulomb interaction, 
Vi ~ 0.5ti, would account for an exchange interaction 
Jeff = 1.2 • 10 3 K which is compatible with the high- 
temperature experimental observation J cxp w 1.4 • 10 3 K. 

IV. SUMMARY 

Using the DMRG method, we provided numerically 
exact results for the spin excitations, the CDW order 
parameter, and the Tomonaga-Luttinger parameter of 
the one-dimensional dimerized extended Hubbard model 
at and near commensurate fillings. 

In the presence of a dimcrization we confirm numer- 
ically that gap for the spin excitation is finite at half 
band-filling. However, the gap immediately disappears 
when the system is doped infinitesimally because there is 
no mechanism which confines the holes to a single dimer. 
This result is qualitatively consistent with a rapid sup- 
pression of the spin gap with Zn doping in the spin-Peierls 
Heisenberg system CuGeOg^,, irrespective of the diffi- 
culty in metalization this material^. 

For the Tomonaga-Luttinger parameter the effects of 
the dimerization are weak in the absence of the nearest- 
neighbor Coulomb interaction V and away from quarter 
filling. At and near quarter filling, the lower Peierls band 
is essentially half filled and the dimerized Hubbard model 
at filling n = 1/2 ±5 can be understood qualitatively and 



even semi-quantitatively in terms of an effective single- 
band Hubbard model at electron density 2n. From the 
result of the corresponding Hubbard model at half band- 
filling it immediately follows that K p = 1/2 holds for 
the dimerized Hubbard model at infinitesimally doping 
away from quarter filling. Therefore, the Tomonaga- 
Luttinger parameter for the weakly doped quarter-filled 
system sensitively depends on the strength of the dimer- 
ization. In general, the dimerization tends to reduce 
K p gradually because the effective scattering processes 
within the Peierls bands increase with the size of the 
Peierls gap. 

In the presence of the nearest-neighbor Coulomb in- 
teraction, the case of quarter filling also deserves special 
attention because the Mott-Hubbard insulator goes over 
to a CDW insulator with a finite spin gap at a critical 
interaction strength V c . The dimerization opposes the 
formation of the CDW phase, for example, the critical 
nearest-neighbor interaction shifts from V c /t\ « 2.65 in 
the absence of dimerization to V c /t\ = 4 in the dimer 
limit. 

The suppression of the charge order at quarter filling 
by the dimerization is reflected in a tendency to stabilize 
the metallic state by the dimerization away from quar- 
ter band-filling. However, the increase of the electron- 
electron scattering by the nearest-neighbor Coulomb in- 
teraction overcomes that tendency and results in a net 
reduction of K p as a function of the dimerization and 
the nearest-neighbor interaction, see Fig. [7| As a con- 
sequence, fairly small values, K p w 0.25, can be ob- 
tained for a moderate five-percent doping of the quarter- 
filled dimerized extended Hubbard model at moderate 
Coulomb couplings, U/ti — 6, Vjt\ = 3. 

It is difficult to reconcile all experimental data for the 
Bechgaard salts with our findings for the dimerized ex- 
tended Hubbard model in one dimension. In order to find 
small values for the Tomonaga-Luttinger parameter, the 
Coulomb interactions must be large enough to reach the 
region of a (doped) CDW insulator which is not easily 
reconciled with the high-temperature data for the ex- 
change interaction. We suspect that the one-dimensional 
dimerized extended Hubbard model is still too simplistic 
to describe the physics of the Bechgaard salts adequately. 
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